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An excess in 7 -rays from the galactic center observed by the FERMI telescope has been proposed 
as a possible signal of dark matter annihilation. Recently, the FERMI collaboration showed that sys¬ 
tematic errors broaden the range of spectral shapes for this excess. We demonstrate fits to this range 
for ( 1 ) flavor-violating annihilations to top-charm pairs and ( 2 ) annihilations to on-shell bosonic 
mediators which decay to Standard Model quarks in a boosted frame. Annihilation of 40 - 100 GeV 
DM to pairs of spin-1 mediators provide a good fit to the FERMI spectrum with a normalization 
consistent with a thermal relic. Top-charm modes and annihilation to three pseudoscalar mediators 
can fit the spectral shape but typically require non-thermal annihilation cross sections. 


I. INTRODUCTION 

Analyses of the data from the FERMI Space Telescope 
show an excess of 7 -rays with energy 1-10 GeV ema¬ 
nating from the center of the galaxy [ 1 - 11 ] and at high 
galactic latitudes [12-14]. More recent studies, including 
one by the FERMI collaboration [15], include estimates of 
the systematic uncertainties on this excess [ 10 , 11 , 16] 
and have broadened the range of allowed spectra. Pro¬ 
posed astrophysical explanations exist [17-21], but it is 
not yet clear if any single astrophysical mechanism can 
account for the excess on all angular scales. 

An intriguing alternate origin for the excess is the anni¬ 
hilation of dark matter (DM) into Standard Model (SM) 
final states which later shower to photons. Early DM fits 
to the excess preferred 0(40) GeV DM annihilating into 
b pairs or 0(10) GeV DM annihilating into lepton pairs 
with a cross section compatible with that required for a 
thermal relic; see [22-25] for phenomenological models. 
Including the systematic uncertainties from more recent 
studies [15, 16] allows a range of 0(10 — 100) GeV DM 
masses and several options for SM final state pairs [26] (a 
similar range of masses was noted in [4]) and can alleviate 
bounds from dwarf spheroidals [27]. 

When DM annihilates into two SM particles, the 7 - 
ray spectrum is purely determined by the identity of the 
SM final state and the dark matter mass. If these in¬ 
teractions are mediated by heavy particles, the couplings 
required to reproduce the 7 -ray excess are constrained 
by monojet [24] and direct detection experiments (see 
e.g. [28]). This suggests a dark sector that contains light 
mediators. In the limit where these mediators are lighter 
than the dark matter, DM can annihilate to on-shell me¬ 
diators, and the annihilation rate is independent of the 
mediator coupling to the SM [29]. In this way, one may 
parametrically avoid the simplest direct detection and 
collider bounds while preserving the FERMI signal. This 
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FIG. 1. Benchmark 7 spectra for s-wave Dirac DM annihila¬ 
tion at the thermal relic cross section: bb (blue), on-shell vec¬ 
tors going to light quarks (red), on-shell pseudoscalars going 
to bb (green). The yellow band is the envelope in (3). Dashed 
lines correspond to the expected spectral shape without the 
boost from on-shell mediators (see text). 


scenario was explored recently in [30-33] for annihilation 
into quarks, [34] for annihilation into electrons, and [35] 
for the case of multistep cascades introduced in [36]. 

In this note, we compare the excess presented by the 
FERMI collaboration [15] to three new types of dark mat¬ 
ter models. In the first, dark matter annihilates to a 
tc final state. In the second, dark matter annihilates to 
two vector mediators which each decay to two quarks. In 
the third, dark matter annihilates to three pseudoscalars, 
which each decay to two quarks. For each class of mod¬ 
els, we use the pppc package [37, 38] for Mathematica [? 
] to generate spectra. We compare these to the spec¬ 
tra presented by the FERMI collaboration, and find the 
allowed parameter space of the model. 

We find that many of these models are viable; for in¬ 
stance, as shown in Fig. 1, the spectra are well fit by the 
decay to vector mediators with a somewhat larger DM 
mass and where the dark matter has the correct thermal 
annihilation cross section. We also show that the shapes 
of the spectra from the two other classes are compatible 
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with the spectra presented by the FERMI collaboration, 
but require a non-thermal component to the relic density. 
Figs. 2 , 5, and 7 summarize our results. We show that 

1. Top-charm annihilation modes provide viable spec¬ 
tra that are very similar to bb. 

2. Vector mediators with universal quark couplings 
are consistent with the spectrum and thermal relic 
abundance for DM masses between 40 - 100 GeV. 

3. Pseudoscalar mediators have a range of masses 
which fit the spectral shape but typically require 
(av) larger than the thermal relic value. 

These qualitative results are only weakly dependent on 
the mediator mass. 


II. 7 SPECTRA FROM DM ANNIHILATION 

The spectrum of the photon flux at Earth d^/dE from 
the annihilation of DM y of mass and with galactic 
density profile Pdm(^) is 

r r^drdQ pI^ {av) dN _ {av) dN ^ 
dE J 47 rr^ r] dE 47rrjm‘^ dE 

where (av) is the thermally averaged cross section for 
DM annihilation and dN/dE is the spectrum of photons 
from the final states and is encoded in pppc [37]. The 
integral of the squared DM density profile is J = 1.08 x 
10^^ GeV^/cm^ assuming an NFW profile [39, 40] with 
7 = 1.2 following the 15° x 15° region of interest about 
the galactic center and fit by the FERMI analysis [15]. r] = 
2(4) for Majorana (Dirac) dark matter^. One typically 
assumes that annihilation is s-wave and consistent with 
the rate required for a thermal relic [42, 43], 

(av) ^ l.lf] X 10“^^ cm^/s. (2) 

For the remainder of this paper we assume Dirac DM so 
that ?7 = 4. In the case of on-shell mediators the spectra 
dN/dE are (1) smeared out due to the boosting of the 
SM states produced from mediator decay ( 2 ) enhanced 
because each annihilation produces a multiplicity of final 
states according to the number of mediators; see [32] for 
details. 

We compare the DM annihilation spectra of on-shell 
mediator models to those observed by the FERMI tele¬ 
scope. We refer to [26] for a fit to models with DM an¬ 
nihilation to pairs of SM final states. However, the anal¬ 
yses in [ 11 , 16] notwithstanding, it is difficult to quantify 


^ These values have different origins. For Majorana DM, this cor¬ 
rects a double counting in the thermally averaged initial state 
phase space in {av). For Dirac DM, this accounts for annihila¬ 
tion being proportional to See, e.g. [41]. 
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FIG. 2. Comparison of the tc (red) and bb (green) spectra for 
a choice of parameters. The bands show normalizations that 
fit the envelope (3), corresponding cross sections are given: 
the upper (lower) range corresponds to the lighter (heavier) 
DM, shown in green (red). Dashed lines show the shape of 
nearby normalizations. 


systematic errors in this excess. We thus take a comple¬ 
mentary approach and compare our spectra to the power 
law with exponential cutoff fits from FERMI collabora¬ 
tion analysis presented in [15]. We use the four spectra 
E^d^i/dE in that analysis to define an envelope that is 
meant to estimate the systematic error on the excess: 


d^miniE) 

dE 

dE 


= a ^miiij 


= a max.- 


dE J ’ 


(3) 


where a a scaling of the envelope size that we choose to 
be 2 in this paper^; see Fig. 1. We quantify the fit of DM 
annihilation spectra relative to the envelope (3) in two 
ways: 


1. For a fixed envelope normalization a, we specify 
the range of annihilation cross sections {av) that 
fall within the envelope. This makes it simple to 
compare to the target thermal relic cross section. 

2. For a given final state spectrum dN/dE, we specify 
the minimum envelope normalization a that may 
accommodate this shape. This quantifies the qual¬ 
ity of the spectral shape fit. We present these plots 
in Appendix A. 


We emphasize that this is an estimate of the parameter 
space for the on-shell mediator scenario that is comple¬ 
mentary to the conventional y^ analysis in [26]. 


^ This is comparable to the range J' G [0.14,4.0] in §3.2 of [26]. 
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FIG. 3. Annihilation cross section required for tc mode to fit 
within the envelope (3). 


III. ANNIHILATION TO TOP-CHARM 
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Previous analyses of DM annihilation to pairs of SM 
particles focus on flavor-conserving fermion pairs or bo¬ 
son [44]. One may also consider flavor-violating modes. 
For final state particles i and j with mass splitting, 
Ami = ml — the injection energies differ, 


Ei = 


4m^ -h 
4m^ 


( 4 ) 


The resulting photon spectrum is a combination of two 
spectra with different characteristic energy scales. The 
mass differences for first and second generation quarks 
are too small for this to be an appreciable effect. Fur¬ 
ther, transitions between 6-quarks and lighter quarks are 
tightly constrained from the 5-factories [45, 46]. Thus 
the most plausible charge-neutral combination is a tc in¬ 
teraction [47, 48]; this has a range of UV motivations 
and collider signatures that are unique from the flavor- 
conserving case [49-51]. We leave lepton flavor-violating 
modes for future work [52]. 

Note that the ti spectrum is generally a difficult fit 
to the FERMI excess due to the hardness of the spec¬ 
trum [26]; tc alleviates this since it only requires 2mx ^ 
180 GeV rather than 2mt and this can better fit the enve¬ 
lope (3). We show this by comparing to a yy ^ bb spec¬ 
trum in Fig. 2. Note the similarity of the two spectra for 
the choice of parameters. Fig. 3 shows the annihilation 
cross section required to fit in the envelope with a = 2; 
one is restricted to values of near its lowest allowed 
value to be near the thermal cross section. 


IV. ANNIHILATION TO SPIN-1 MEDIATORS 

The leading 5-wave contribution to DM annihilation to 
on-shell spin-1 mediators is yy ^ VV, where V may be 
either CP even or odd [30, 32, 33]. With the modest as¬ 
sumption of minimal flavor violation, one expects such a 
mediator to couple universally to all quark flavors. Here 


FIG. 4. Separation of the (av) required to fit in the FERMI en¬ 
velope from the thermal value in units of the thermal value for 
the case of annihilation to two spin-1 mediators which decay 
to light quarks. The range G [50,100] GeV is consistent 
with a thermal relic. 
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FIG. 5. Extreme spectra for annihilation to on-shell, univer¬ 
sally coupled vector mediators. Gonventions as in Fig. 2. 


we present results for universal couplings to light quarks; 
the modification from the b quark contribution is a per¬ 
cent level correction. For completeness we present results 
a 6-philic vector in Appendix C. 

As noted above, the on-shell mediator smears out the 
photon spectrum. In the non-relativistic limit, kinemat¬ 
ics restrict the Vs to be mono-energetic. However, the 
photon spectra from V ^ qq must be convoluted with a 
box spectrum that encodes the boost from the mediator 
rest frame [32]. This smears out the spectrum, as demon¬ 
strated by the dashed lines in Fig. 1. Fig. 5 plots the 
shapes of two extreme examples of allowed DM masses 
in this scenario. The range of normalizations that allow 
each shape to fit the envelope (3) are listed for compari¬ 
son with (2). 

This scenario is able to accommodate a thermal relic 
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FIG. 6. Separation of the {av) required to fit in the FERMI 
envelope from the thermal value in units of the thermal value 
for the case of annihilation to three pseudoscalar mediators 
which decay to bs. The thermal cross section is generally not 
large enough not explain the FERMI excess. 


within the FERMI envelope (3). We note that this class 
of mediators is constrained theoretically by anomaly can¬ 
cellation [53]. See [54] for a recent exploration of viable 
Z' models. 



Photon Energy [GeV] 


FIG. 7. Extreme spectra for annihilation to on-shell, 5-philic 
pseudoscalar mediators. Gonventions as in Fig. 2. 



DM mass, [GeV] 


V. THREE PSEUDOSCALAR MEDIATORS 


FIG. 8. Minimum a required in (3) for yy ^ tc. The spectral 
shape is in tension with the envelope for ^180 GeV. 


For spin-0 mediators, the leading 5-wave on-shell an¬ 
nihilation mode comes from the annihilation to three 
scalars [32]. We consider the simplest case, xx ^ for 
a pseudoscalar (p. The energy spectrum of these media¬ 
tors is no longer a simple box-spectrum and must be con¬ 
voluted with the photon spectra. We present expressions 
for calculating this mediator spectrum in Appendix B. In 
contrast to spin-1 mediators, minimal flavor violation im¬ 
poses that spin-0 mediators couple proportionally to the 
Yukawas so we make the approximation of an exclusive 
if ^ hh decay. 

We show fits for this scenario in Figs. 6 and 7. There 
is some tension with the thermal relic cross section. This 
can be seen in the expression for the photon flux (1): 
while dN/dE is three times larger than that of XX ^ bb 
due to the multiplicity of 6s, the flux also decreases by 
9 due to the more sparse DM number density associated 
with a heavier x- This is because must be approx¬ 
imately three times heavier so that each b has roughly 
the same injection energy as in the xX ^ bb case. Hence 
{av) must compensate for this overall 0(1/3) factor. In 
fact, as explained in [17], this is still an over-estimate 
for comparison to the thermal relic cross section since at 
freeze out, the p-wave annihilation mode to two cps is not 
significantly suppressed relative to the 30 mode. Thus, 


this case may require a non-thermal DM abundance, see 
e.g. [55-57]. 

A template simplified model for pseudoscalar media¬ 
tors and the FERMI excess was presented in [58] with the 
on-shell limit considered in [32]. Since spin-0 particles 
couple left- and right-chiral fermions, UV completions 
of this framework typically require interactions with the 
Higgs [59]; see recent studies in two-Higgs doublet [60] 
and (N)MSSM [61-63] frameworks. 
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Appendix A: Shape Fits 


an energy function 


As explained in Sec. II, there are two ways to quantify 
the fit of a particular photon spectrum. In the main text, 
Figs. 4 and 6 showed contours of {av) required for a given 
spectrum to fit in the envelope (3). This quantifies overall 
normalization. Alternately, one may quantify the fit of 
the spectral shape by the minimum envelope scaling a 
for which there exists any {(jv) that allows the spectrum 
to fit the envelope. This data is plotted for the cases in 
the main text in Figs. 8 and 9. 


Appendix B: s-wave Pseudoscalar Spectrum 


We summarize the spectrum dN^jdE^ of non¬ 
relativist ic DM annihilation into three pseudoscalars. For 
simplicity we set the coupling to one. The squared am¬ 
plitude is 

4 1*^1 ~ ^^^bRcd) A^B.C^D ^ 

spins ABCD 


where A^B^C^D G {1, 2, 3} label the final state if and 


Rab 


{P% + m^) {-Qa + m^) + Pa • Qb 

(m2 - 2mEA) (m^ - 2mEB) 


(B2) 


written with respect to the 4-vectors formed from the 
outgoing f momenta Pa and the incoming DM (anti-DM) 
momenta kiff 


II 

1 

to 

1 

II 

0 

(B3) 

The differential cross section is 


II 

to 

1 pidpi p2dp2 

(B4) 

(27r)3 2Ei 2 E 2 ■ 

In order to find d{av) jdE^ 

= ‘^d{(jv)/dEi^ we must in- 


tegrate over dp 2 over its allowed kinematic range. Define 


e(Pi,P2,cos0) = El+E 2 + E3{pi,p2,cos0), (B5) 

where 9 is the angle between pi and p 2 , pi ,2 are lengths 
of 3-momenta, E 12 = Pi ,2 + ^^5 

^3 = Pi + P 2 + 2pip2 cos 0 (B6) 

The region of dp2 integration is given by the con¬ 
dition that there exists some |cos6>| < 1 such that 
e(pi,p 2 , cos = ^/s = 2m^. This can be visualized by 
plotting e(pi,p 2 5 — 1) as a function of p 2 and noting that 
for general values of cos6>, e lies above this line. There 
are two regions. First, when e(pi, 0, — 1) > then pf^"^ 
and pf^^ are given by the solution to 

= (B7) 

The second region has e(pi, 0, cos6>) < y^, in this case 
pf^"^ and pf^^ are given respectively by the solutions to 

e(pi,p™“,+l) = x/^ (B8) 

= (B9) 

Upon performing the dp2 integration, the f spectrum can 

be found using the density of states [64], 

J_^ ^ 1 d{av) 

% dE^ {av) dE^' ^ ’ 


Appendix C: MFV-Violating Modes 

In the main text we made the point that minimal fla¬ 
vor violation (MFV) suggests that spin-1 mediators cou¬ 
ple flavor-universally while spin-0 mediators couple ac¬ 
cording to the Yukawas. For completeness, we present 
summary plots of the opposite scenarios in Fig. 10. 
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FIG. 9. Contour plots for the minimum envelope scaling to accommodate the spectral shape of on-shell spin -1 (left) and 
pseudoscalar (right) mediators. Contours correspond to a in (3). 
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FIG. 10. Summary plots for MFV-violating couplings following the figures in the main text. 
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